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ABSTRACT 

We investigate torsional Alfven modes of relativistic stars with a global dipole mag- 
netic field. It has been noted recently (Glampedakis et al. 2006) that such oscillation 
modes could serve as as an alternative explanation (in contrast to torsional crustal 
modes) for the SGR phenomenon, if the magnetic field is not confined to the crust. 
We compute global Alfven modes for a representative sample of equations of state 
and magnetar masses, in the ideal MHD approximation and ignoring ^ ±2 terms in 
the eigcnfunction. We find that the presence of a realistic crust has a negligible effect 
on Alfven modes for B > 4 x 10^^ G. Furthermore, we find strong avoided crossings 
between torsional Alfven modes and torsional crust modes. For magnetar-like mag- 
netic field strengths, the spacing between consecutive Alfven modes is of the same 
order as the gap of avoided crossings. As a result, it is not possible to identify modes 
of predominantly crustal character and all oscillations are predominantly Alfven-like. 
Interestingly, we find excellent agreement between our computed frequencies and ob- 
served frequencies in two SGRs, for a maximum magnetic field strength in the range 
of (0.8-1.2) xlQif' G. 
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1 INTRODUCTION 

Soft Gamma Repeaters (SGRs) have attracted a lot of attention recently, since there exist at least two sources in which 
quasi-periodic oscillations have been observed. The frequency of these oscillat ions, in the tail of the burst , indicate that they 
must originate in the vicinity of the SGR source, thought to be a magnetar ( Duncan &: ThompsonI 1992) . A popular m odel 
for explaining these oscillations involves the excitation of torsional mo des in the solid crust of the star ( Duncan 19981 ) . An 
alternative explanation has been provided by Glampedakis et al. ( 20061 ) . In a simple toy model, they showed that, in addition 
to crustal modes, a magnetized star could also have a discrete spectrum of global Alfven modes, which would dominate over 
crustal modes for sufficiently strong magnetic fields. Here we investigate global, torsional Alfven modes in realistic models of 
relativistic stars with a global dipole magnetic field and compare our computed frequencies with observations. 

SGRs produce giant flares with peak luminosities of 10** - 10*® erg/s, which display a decaying tail for several hundred 
seconds. Up to now, three giant flares have been detected, SGR 0526-66 in 1979, SGR 1900-1-14 in 1998, and SGR 1086-20 
in 2004. The timing analysis of the latter two events revealed several quasi-periodic oscillations (QPOs) in the decaying tail, 
with frequencies of 18 , 26, 2 9, 92.5, 150, 626.5, and 1837 Hz for SGR 1806-20 and 28, 54, 84, and 155 Hz for SGR 1900-H4 (see 
Strohmaver fc Watts! (|2006l ) and references therein). Furthe rmore, for SGR 1806-20 the possible detection of two additional 
higher frequencies (720 Hz and 2384 Hz) has been reported ( Strohmaver fc Watts! iooil ) . 

Spherical stars have generally two type oscillations, spheroidal modes with even parity and the toroidal modes with axial 
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parity. The observed QPOs in SGR tails are thought to originate from toroidal oscillations, since the latter are excited more 
easily than poloidal oscillations, as toroidal modes do not involve densit y variations. In Newton ia n theory, there ha v e bee n 
several invest i gation s of torsional mode s inneutron star crusts, see e.g. Hansen fc CiofB ( 1980l ): McDermott et al. ( 19881 ): 

O n the other hand, o n ly a few studies ha v e taken general relativity 



Carroll et al 



into account 



( 1986h: Istrohmaveil (Il99lh : 



Schumaker fc Thornel ijlQSi ): 




.ed (12001 
LeinsI (1 199 



Messios et al.l 



(200a). The more recent of these studies (see e.g. ISotani et al.l ( 20061 ) . hereafter Paper I) have shown that some of the 



200ll ): ISotani et all (|2006h : ISamuelsson fc AnderssorJ 



observational data of SGRs could agree with the crust torsional modes, if, e.g., frequencies lower than 155 Hz are identified 
with the fundamental modes of different harmonic index £, while higher frequencies are identified with overtones. However, as 
mentioned in Paper I it will be quite challenging to identify all observed QPO frequencies with crustal torsional modes. For 
example, it is difficult to explain both the frequencies of 26 and 29 Hz for SGR 1806-20 for a single stellar model, because the 
actual spacing of torsional modes of the crust is larger than the difference between these two frequencies. Similarly, the spacing 
between the 626.5 and 720 Hz QPOs in SGR 1806-20 may be too small to be explained by consecutive overtones of crustal 
torsional modes. In the present paper, we show that such difficulties are avoided if one identifies the observed frequencies with 
global Alfven modes. 

The article is structured as follows: in the next section [2] we briefly summarize the equations to be solved and the 
assumptions made. In section |3] we discuss the numerical results, where to understand qualitatively the properties of global 
Alfven modes we also treat two special cases. The comparison with the observed data in SGRs are done in section |4] and 
finally we give a summary and discussion in section [5] Unless otherwise noted, we adopt units of c = G = 1, where c and G 
denote the speed of light and the gravitational constant, respectively, while the metric signature is (— , +, +, +). 



2 MAIN EQUATIONS AND NUMERICAL METHOD 



The main equations and numerical method required for the present study have already been derived in Paper I. Here, we only 
give a brief summary of the main assumptions and final relations. 

We assume that the equilibrium stellar model is described by a solution of the TOV equations for nonrotating, relativistic 
stars and a metric of the form 



ds^ — — e'^*df^ + e^^dr^ -f r^{d6^ + sin'^ 9d(j)^), 



(1) 



where $ and A are functions of the Schwarzschild radial coordinate r. We neglect the influence of the magnetic field on 
the equilibrium configuration, since the magnetic energy £^4 is negligible, compared to the gravitational binding energy £g, 
Sm/Ss « 10-\B/{W^<'G)f. 

Adopting the ideal MHD approximation, we assume that the star is endowed with a dipolar magnetic field, described by 



He = 



-ai . 



47r 



-ai, 



(2) 
(3) 



where //^ = B^^jy^, with B^^ the magnetic field 4- vector, and a\ being a function describing the radial behavior of the 
magnetic field. In the exterior of the star 



(cx) 



'8M3' 



In 



/ 2M\ 2M 
\ r ) r 



2M 2M^ 



+ 



(4) 



where fj,b is the magnetic dipole moment for an observer at infinity ( Wasserman fc Shapiro 19831 ). while the interior solution 
is found by solving Maxw ell's equations for the 0-component of the electromagnetic 4-potential A,^, assuming a particular 



form for the 4-current Jj, (Konnoetal 



19991 ). 



We have used a variety of neutr on star models, using fo ur different equations of state fo r the core, ranging from a very soft 
EoS (E oS A) JPandharipandj 197ll ) to a very stiff (EoS L) (jPandharipande fc Smithl 19751 ). with two intermediate ones, EoS 
WFF3 ( Wiringa et al.l 1988 ) and APR ( Akmal et al. 19981 ). For each EoS we constructed a number of models, starting from 
a gravitational mass of 1.4M0 and reaching close to the maximum mass limit in increments of O.2M0. In order to separately 
investigate the effect of the composition of the crust, we matched the various high-density EoS to two different proposed 
equations of state for the crust, one recent derived bv lDouchin fc Haensell (|2001f ) (DH) and, for reference, and older EoS by 
Negele fc VautherinI (119731 ) (NV). The two crust EoS differ signific antly both in the de t ailed composition, as well as in the 
de nsity at the base of the crus t, which is at p ~ 2.4 x lO^'^gr/cm"' for Negele fc VautherinI (1973) and at p ~ 1.28 x lO^^'gr/cm^ 
for Douchin fc Haensell (|200ll ). For the stiff EoS these different properties of the crust have a considerable effect on the bulk 
properties. Figure 1 of Paper I displays the mass-radius relationship of nonrotating equilibrium models constructed with the 
various high-density EoS, in combination with our two choices for the crust EoS. The individual models for which we compute 
torsional modes are shown with symbols, while their detailed properties are also listed in Table 1 of Paper I. 
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For studying individual oscillation modes, the linearized dynamical equations governing the magnetized fluid and crust 
are specialized to the case of axial perturbations, adopting the relativistic Cowling approximation (i.e. neglecting spacetime 
perturbations). In this case, the only nonvanishing perturbed fluid variable is (assuming a harmonic time dependence with 
frequency ui) 

Su'*' = iLje--^y{r)e''^%{e), (5) 
where 

1 



m 



-dePeicose), 



(6) 



and de denotes the partial derivative with respect to 9. Above, y{r) describes the radial dependence of the angular displacement 
of the stellar material, Pe{cosO) is the Legendre polynomial of order i and we have set m. = 0, due to the degeneracy in m for 
a spherically symmetric background. 

The linearized induction equation in the MHD approximation relates the perturbation in the magnetic field to y{r). 
Neglecting magnetic-field-induced i±2 couplings leads to the eigenvalue equation 



^ + (H-2Ai) 



ai 



y" + {(^ + $'-A')M + M' + (l + 2Ai)-^[(c&'-A')ai-f2ai']}y 



+ 



e + (1 + 2Ai) 



ai 



\ '2 



-(A -2) 



lie 



27rr4 



+ (2 + 5A02^{(<l>'-A')a.'-fa/'} = 0, 



(7) 



where fi is the shear modulus in the crust and 
A = £{£+1), 
Ai = 



£{£+!) 



{2£-l){2£ + 3)' 

In order to solve a system of first-order ODEs, we define new variables and 3^2, through 

yi = yr'-', 
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(l + 2Ai 
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2 2{A-*) 

uj re • ' 
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lie 



2TTr^ 



+ (2 + 5Ai) 



aie 



e y\ 3^2, 



(8) 
(9) 

(10) 

(11) 



In defining the new variables, we took into account the form of the eigenfunction y near the center, where y ^ ^ . The 
final first-order system of equations to be solved numerically for the real eigenvalues ui is then 



(12) 



(13) 



where the unperturbed Maxwell's equations were used in order to eliminate the term of a\" . 

The above system of equations is solved as an eigenvalue problem, imposing a) regularity at the center, b) a continuous 
traction condition at the crust-core interface, and c) a zero-torque condition at the stellar surface. Regularity at the center 
implies 

^ (14) 



3^2 = (^ - 1) 



^i + (l + 2Ai) — 



e 3^1 . 



The demand that the traction is continuous at the crust-core interface implies 



y 



'(-) 



1 + 



y 



'(+) 



(15) 



where Us and ua are the shear and Alfven velocities, respectively. With our choice of variables, th is condition becomes 
y^^^ = 3^2^'' . At the stellar surface, the zero-torque condition ST^"^"^ — ( Schumaker fc Thornel 1983 ) implies 3^2 = 0. We 
neglect the possible presence of a thin fiuid ocean. 
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Figure 1. Frequencies of Alfven modes, fa„, with 1 = 2 and n = 0, 1, and 2 as a function of the normaUzed magnetic field B/B^. The 
equilibrium models have mass 1.4Mq and only results for a very soft (A) and a very stiff (L) high-density EOS are shown, combined with 
the DH EOS for the crust (see text). Individual numerical results are shown with various marks, while the continuous lines correspond 
to the empirical formula l|17| l with coefficient values in shown in Tables [T] [2] and|3l 

Specific frequencies of predominantly or pure crustal modes are labeled as itn, while the frequencies of predominantly or 
pure Alven modes are labeled as tan, where I is the angular index and n is associated with the number of radial nodes in the 
eigenfunction of various overtones. 



3 NUMERICAL RESULTS 

Our equilibrium models comprise both a fluid part and a thin solid crust, which can have a thickness ranging from about 3% 
to 12% of the stellar radius, depending on the EOS and mass. We will compute global torsional Alfven modes in the presence 
of both of these components. However, in order to understand qualitatively the properties of these modes, we will first study 
them in the limiting case where one neglects the presence of a solid crust and treats the model as a pure magneto-fluid. 
In addition, in order to understand the influence of the presence of a solid crust on the Alfven modes (and, conversely, the 
influence of a global magneto-fluid on the torsional crust modes) we will study a toy model of varying crust thickness, before 
focusing on the realistic models presented in the previous section. 

3.1 Global Alfven modes for a pure magneto-fluid 

As a limiting case, we flrst consider global Alfven modes, in the absence of a solid crust, i.e. setting = in the crust region. 
These modes are pure Alfven modes and degenerate to zero frequency in the limit of a vanishing magnetic field, since then 
the eigenvalue problem becomes 



which has = as a trivial solution. In the presence of a magnetic field, the eigenfrequency becomes proportional to the 
magnetic field strength B. 

Figure [T] shows the frequency of several modes (with 1 = 2 and n = 0, 1, 2) as a function of the strength of the magnetic 
field, for IAMq models of a very soft (A) and a very stiff (L) EOS. Various lines in the figure are linear fits to the numerical 
data. As is apparent from these numerical results, the linearity of the frequency as a function of B is preserved to a very high 
accuracy, at least in the range of B = 10" - 10" G. Therefore, our numerical results can be represented by the following 
empirical formula 



(e + p)e 



-2*+2A 2 



i^"y = 0, 



(16) 




(17) 
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Table 1. The values for the fitting factors if^o (Hz) of equation 
I I17I I. The fitting factors have been calculated for magnetic field 
strengths up to 10^^ G. 
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where is a typical magnetic field strength, which we take to be Bp = 4 x 10^^ G (as in Paper I) and i(3n are coefficients, 
which depend on the equilibrium model and the mode in question. For each equilibrium stellar model the coefficients ePn for 
^ = 2 - 8 and n = 0, 1, 2 are shown in Tables [T] [21 and El 

As a general trend, one sees that the value of the coefficients il3„ becomes larger (and thus the dependence of the mode 
frequency on the magnetic field strength becomes stronger), as the mass of the star increases (for a given EOS). The same 
trend holds with increasing stiffness of the EOS, for a given mass. One also sees from the Table[T]that f/3o/2/3o ~ \/ so 
that the following empirical formula holds for the fundamental modes 

Examining more closely the dependence of the mode frequency on the parameters of the equilibrium models, we find that 
the coefficients f/3o are a linear function of the inverse of the compactness M/R (see Figure [2]|, so that the frequency of the 
fundamental modes is approximately given by 




(with an accuracy of a few per cent). For given strength of the magnetic field the frequency of the fundamental modes varies 
by a factor of up to 2.4 within the set equilibrium models in Table [T] Similarly, the frequency of the n = 1 and n = 2 
overtones of the 1 = 2 mode varies by up to a factor of 2.7. It follows that the frequency of global Alfven modes is sensitive to 
the magnetic field strength and the bulk properties of the equilibrium model and a successful identification of several modes 
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Table 2. The values for the fitting factors if^i (Hz) of equation 
I I17I I . The fitting factors have been calculated for magnetic field 
strengths up to 10^^ G. 



Model 


2/3l 


3/3l 


4/3l 


5/3l 


6/3l 


7/3l 


8/3l 




'^7 09 
1 .uz 


42 91 


A7 9fi 


QA 


54 24 


^7 9fi 


fin 1 




ZD. ( i 


Ol .UD 


04t.0U 


•^7 fiH 


oy .oy 


41 09 

^1 .yz 


AA 1 A 


WFF3+DHi4 


44.50 


51.45 


56.54 


60.83 


64.65 


68.17 


71.48 


WFF3+DHi6 


35.80 


41.48 


45.67 


49.22 


52.39 


55.32 


58.08 


WFF3+DHi8 


26.46 


30.79 


34.04 


36.81 


39.32 


41.65 


43.87 


APR+DH14 


51.63 


59.53 


65.28 


70.12 


74.42 


78.39 


82.11 


APR+DH16 


43.51 


50.27 


55.21 


59.37 


63.09 


66.51 


69.73 


APR+DH18 


36.56 


42.32 


46.55 


50.13 


53.34 


56.30 


59.09 


APR+DH20 


30.30 


35.15 


38.74 


41.79 


44.53 


47.07 


49.47 


APR+DH22 


24.18 


28.14 


30.43 


33.63 


35.92 


38.05 


40.08 


L+DH14 


61.38 


70.83 


77.71 


83.48 


88.60 


93.29 


97.66 


L+DH16 


53.30 


61.56 


67.57 


72.61 


77.08 


81.17 


84.99 


L+DH18 


46.51 


53.78 


59.08 


63.53 


67.47 


71.08 


74.45 


L+DH20 


40.55 


46.93 


51.61 


55.54 


59.03 


62.23 


65.22 


L+DH22 


35.18 


40.74 


44.84 


48.29 


51.36 


54.19 


56.84 


L+DH24 


30.13 


34.93 


38.47 


41.46 


44.14 


46.61 


48.94 


L+DH26 


24.75 


28.77 


31.75 


34.29 


36.57 


38.69 


40.70 


A J_1S.T'\/'-. ^ 
V 14 


00.00 


AO 7^^ 


A7 (17 


f^n 7A 


^A 09 


^7 flf^ 


f^Q on 
oy .yvj 


A -UMV-, ^ 


26 62 


'in 

ou. yo 


'^A ^Q 




'^Q A^ 

oy 


A^ 77 


A'^ QS 


WFF3+NV14 


44.43 


51.37 


56.45 


60.73 


64.55 


68.06 


71.36 


WFF3+NV16 


35.77 


41.44 


45.63 


49.16 


52.33 


55.26 


58.02 


WFF3+NV18 


26.43 


30.76 


34.00 


36.77 


39.27 


41.60 


43.82 


APR+NV14 


49.79 


57.43 


63.00 


67.66 


71.81 


75.61 


79.18 


APR+NV16 


42.12 


48.66 


53.44 


57.46 


61.04 


64.33 


67.42 


APR+NV18 


35.53 


41.13 


45.23 


48.70 


51.79 


54.64 


57.32 


APR+NV20 


29.56 


34.29 


37.78 


40.74 


43.39 


45.85 


48.16 


APR+NV22 


23.68 


27.54 


30.43 


32.90 


35.13 


37.20 


39.16 


L+NV14 


56.55 


65.42 


71.91 


77.35 


82.18 


86.59 


90.69 


L+NViB 


49.07 


56.88 


62.60 


67.40 


71.65 


75.52 


79.12 


L+NV18 


42.73 


49.61 


54.67 


58.92 


62.69 


66.13 


69.32 


L+NV20 


37.20 


43.21 


47.65 


51.40 


54.73 


57.77 


60.60 


L+NV22 


32.28 


37.50 


41.36 


44.62 


47.52 


50.18 


52.67 


L+NV24 


27.68 


32.15 


35.45 


38.25 


40.75 


43.05 


45.20 


L+NV26 


22.68 


26.39 


29.15 


31.49 


33.58 


35.52 


37.34 



would, in principle, allow the inference of both the magnetic field and the compactness of the star, through empirical relations 
such as Equation (|19|l . 

Figure [3] shows the frequencies with £ = 2, 3, and 4, and n = 0, 1, and 2 as a function of the magnetic field strength 
B/Bn, for the particular stellar model A+DH14. The Figure demonstrates that the frequencies of overtones depend on the 
value of £, in sharp contrast to the case of pure crust modes, where the frequencies of overtones are almost independent of £ 
(see, e.g. Paper I). 



3.2 Avoided crossings between global Alfven modes and torsional modes of the crust 

How does the presence of a solid crust affect the global Alfven modes? And, conversely, how are the torsional modes of the 
crust affected when global Alfven modes are also present? In order to answer these questions, we construct a toy model in 
which we introduce a crust of varying thickness Ar/R. We focus on an equilibrium models with M = IAMq and R — 10.4 
km, constructed with a relativistic polytropic EOS with polytropic index TV — 1.0 and polytropic constant K = 130 km^. In 
this toy model, the shear modulus is simply taken to he fi — pvg, where = 1-0 x 10* cm/s 

If one assumes that the whole star is solid {Ar/R = 1.0) then the frequency spectrum of torsional modes only consist 
of crust modes. Qualitatively, the dependence of the mode frequencies on the strength of the magnetic field is similar to the 
results in Paper I, i.e. there is a weak dependence up to B/Bf^ ~ 1 (nearly horizontal branch), followed by a rapid increase 
in frequency for larger magnetic field strengths (magnetic-field-dominated branch). When the crust thickness is reduced to 
Ar/R < 1, global Alfven modes are introduced. In such a composite fluid/crust case, the eigenfunctions of torsional crust 
modes penetrate inside the fluid region, i.e. crust modes become global modes and, conversely, the eigenfunctions of Alfven 
modes penetrate the inside crust region. Thus, the allowed discrete global modes are only those for which the eigenfunction 
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Table 3. The values for the fitting factors f /32 (Hz) of equation 
I I17I I. The fitting factors have been calculated for magnetic field 
strengths up to 10^^ G. 
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at the crust/core interface satisfies the continuous traction condition (|15p exactly. As we will see, this strict requirement leads 
to avoided crossings between the two families of modes and to a dramatic change of the mode spectrum in the limit of small 
crust thickness and high magnetic field. 

Figure U shows the first six frequencies for ^ = 2, as a function of the strength of magnetic field, for four different values of 
the crust thickness, Ar/R — 0.95, 0.70. 0.50, and 0.30. As seen Figure|3Ja) even a small fiuid part {Ar/R — 0.95) allows both 
families of modes to be present, e.g. for small magnetic field strengths one can see both horizontal mode sequences (dashed 
lines, corresponding to crust modes) and Alfven-like sequences (continuous lines). However, as the functional dependence of 
the global Alfven modes with increasing magnetic field strength is different than the functional dependence of crust modes, 
the two families of modes appear to cross each other. On closer examination, one sees that, in fact, avoided crossings take 
place. For {Ar/R — 0.95) the gap at avoided crossings is still small, so that both families of modes are clearly visible. In 
addition, since the fiuid region is still small, the Alfven mode frequencies are high and the avoided crossings take place on the 
horizontal branch of the crust modes. 

When the crust thickness is reduced to Ar/R = 0.70, Figure |4{b), the enlargement of the fiuid region leads to a reduction 
of the Alfven mode frequencies. This leads to the avoided crossings taking place at larger values of the magnetic field. The 
avoided crossing of the fundamental mode takes place close to the transition of the crust modes from the horizontal branch 
to the magnetic-field-dominated branch. Thus, a horizontal branch of predominantly crust modes is still present. However, 
the avoided crossings of higher-order overtones take place at the magnetic-field-dominated branch of crust modes. This leads 
to a significant increase in the gap at avoided crossings and one cannot easily distinguish modes that are predominantly crust 
modes. 

Reducing the crust thickness further to Ar/R = 0.50 and 0.30, Figures |4jc),(d), one observes that the avoided crossings 
of low-order modes take place in the magnetic-field-dominated branch of crust modes and the gap in avoided crossings tends 
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Figure 3. The frequencies of Alfven modes for fa„ with £ = 2, 3, and 4 and n = 0, 1, and 2 as functions of the normaUzcd magnetic 
field B/B/i. The solid, dashed, and short-dashed lines correspond to the frequencies of ^ag, ^ai, and ia2, respectively, and for each n the 
frequencies increase with increasing £. For a given £, the frequency of various overtones does not coincide, in sharp contrast with crustal 
torsional modes. The mass of the equilibrium model, constructed with the A+DH EOS, is M = 1.4Mq. 



to become of the same order as the frequency spacing of consecutive Alfven modes. The avoided crossings are now visible 
as only a small local change in the slope of mode-sequences, where one would expect a crust mode to be present. However, 
modes of predominantly crustal character are no longer distinguishable. The same trend continues when the crust thickness 
is reduced even further, so that for realistic values of the crust thickness (3%-12%) the mode spectrum is similar to the pure 
Alfven mode spectrum shown in Figure [S] 

Regarding the eigenfunctions of oscillation modes, at a given magnetic field strength, the lowest-frequency mode corre- 
sponds to a fundamental mode without any node, while the number of nodes is increasing for overtones of higher frequency. 
We find that the position of nodes change at avoided crossings. The inset in Figure |4|b) follows the continuous sequence of 
mode frequencies, which start out as the 211 Alfven mode (first filled circle at B/B^ — 0.1). After an avoided crossing, the 
sequence becomes a fundamental torsional mode, 2to, (second filled circle at B/B^ = 0.2). After a second avoided crossing 
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Figure 4. The behavior of torsional crust modes {it„) and global Alfven modes (;a„) for different toy models, when varying the crust 
thickness Ar/R. The first six eigenfrequencies for 1 = 2 are shown. For comparison, we also plot the case of Ar/R =1.0 (no fluid region) 
as dashed lines in panel (a). 

the continuous sequence changes into fundamental Alfven mode, 200, (third filled circle at B/B^ = 0.3). Figure [5] shows that 
each of the above three eigenfunctions has one node. Specifically, the eigenfunctions for the 20.1 and 2to modes have one node 
in the fiuid region while the 210 mode has one node in the crust. 

We further observe that the amplitude of the 200 and 211 Alfven modes is larger in the fluid region, while the amplitude 
of the 2to mode is larger in the crust region. This is consistent with the interpretation that, even with the presence of 
a large number of avoided crossings, horizontal branches correspond to predominantly crust-like modes, while branches of 
increasing frequency (with increasing magnetic field strength) correspond to predominantly Alfven-like modes (except for the 
magnetically-dominated branch of crust modes). Of course, in the above example a large crust thickness of Ar/R = 0.70 was 
used, where the two family of modes are still distinguishable. For a much smaller crust thickness, the the horizontal branch 
of crust modes cannot be clearly identified (due to the strong avoided crossings), while the magnetically-dominated branch 
of crust modes is merged with the Alfven modes, as both families of modes have nearly the same frequencies and slopes at 
large magnetic field strengths. 

3.3 Global Alfven modes for realistic models 

Realistic stellar models have very thin solid crusts and the frequency spectrum of global modes approaches that of the pure 
Alfven modes. As an example. Figure |6] shows the three lowest-order £ = 2 modes, 210, 2ai, and 212 for model A-I-DH14. 
Pure Alfven modes (no solid crust) are shown as solid lines while the modes in the presence of the crust are shown as dashed 
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Figure 5. Scaled eigenfunctions of the 20.1 (at B/Bf^ = 0.1), 2*0 (a-t B/B^ = 0.2) and 210 (at BjBy, = 0.1) modes for a toy model with 
large crust thickness of Ar/R = 0.70. The number of nodes and the relative amplitude at the crust/fluid regions is consistent with the 
expected behavior for predominantly crustal modes and predominantly Alfven-like modes, respectively. 



lines. The pure crustal fundamental mode (when the magnetic field is confined to the crust region only, see Paper I) is also 
shown as a dash-dotted line. For weak magnetic fields the presence of the solid crust has some effect on the Alfven mode 
frequencies, but, for large magnetic fields the effect diminishes. For larger i, the effect of the presence of a solid crust on the 
mode frequencies persists for somewhat stronger magnetic fields. 

Since the global modes for large magnetic field strengths are almost identical to pure Alfven modes. Equation (|17|l 
describes their frequencies with good accuracy, for large magnetic field strengths. In addition, the non-degeneracy of the 
frequency of overtones for different values of the harmonic index I will be important in trying to match numerical results to 
observed frequencies. 



4 COMPARISON WITH OBSERVED FREQUENCIES IN SGRS 

As mentioned in Paper I, it is difficult to explain all of the observed frequencies in SGRs using pure crustal torsional modes 
(with magnetic fields confined to the crust only). An example is the small separation of the two frequencies of 26 and 29 Hz 
for SGR 1806-20. Here, we find that such frequency data can be easily explained with global Alfven modes, owing to the 
non-degeneracy of overtones for different values of I. 

Figure [7] shows, as an example, various global Alfven mode frequencies for model A-I-DH14. In this figure, the horizontal 
broken lines denote the observational data up to 150 Hz for SGR 1806-20 while the solid lines denote our computed frequencies. 
We find an excellent agreement between observed and computed frequencies for a magnetic field strength of B/B^ ~ 1.25, 
(shown as a vertical dash-dotted line). Specifically, the observed frequencies of of 18, 26, 29, 92.5, and 150 Hz can be identified 
with the 20.0, 400, sflo, 502, and 504 modes, respectively. 

Figure [8] shows a similar excellent agreement of the observed frequencies for SGR 1900-f 14, with the frequencies of the 
200, 700, sfli, and 602, respectively, for the same chosen stellar model A-I-DH14 and a slightly larger magnetic field strength. 
Since the global frequencies are sensitive to both the stellar parameters and the magnetic field strength, an independent 
determination of the compactness of the star would allow for a determination of the magnetic field. 

The higher-frequency observational data for SGR 1806-20, such as 626.5, 720, 1837, and 2384 Hz could be also be identified 
with higher overtones of global Alfven modes, as there exist many higher overtones with different £. For example, for model 
A-I-DH14 and with B/B^^ = 1.25 (as above), we compute the following mode frequencies: 4025 = 632, 4029 — 724, 2<i83 = 1841, 
and 201108 = 23 82 Hz, which are in all good agreement with the observed frequencies. We remark that both the 626.5 and 
720 Hz frequencies can be identified with some global Alfven modes, while this would not be possible with the crustal modes, 
when the magnetic field is confined to the crust (see numerical results in Paper I), due to the near degeneracy of overtones of 
different £. 

Table |4] shows the maximum values of the magnetic field, for which one can at least identify the frequency of the 210 
mode. For models which have no entry in this Table, such an identification was not possible. In reality, some stellar models 
with lower magnetic strength can agree well with the observational data, if the 210 mode is omitted. For example, the stellar 
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Figure 6. The frequencies of global Alfven modes with 1 = 2 and n = 0, 1, and 2 as functions of the normalized magnetic field 
{B/Bfji,). The dashed line corresponds to a realistic model (A+DH14) with both fluid and solid regions. The solid line corresponds the 
same model, but setting = in the crust region (pure magneto-fluid). The differences between the two cases are limited to small 
magnetic field strengths, while for higher values there is coincidence. For comparison, we also show the 2*0 pure crustal mode (when the 
magnetic field is confined to the crust). 




Figure 7. Comparison of several specific Alfven modes (solid lines) with observed frequencies for SGR f 806-20 (horizontal dashed lines) 
for the stellar model A-I-DH14. The agreement for all frequencies up to f50 Hz is excellent, for a magnetic field strength of B/ Bpi 1.25 
(vertical dashed-dotted line). 



model with L-I-NV18 and B/Bf^ = 0.60 can agree well with all observational data lower than 150 Hz, when identified with the 
frequencies of the 400, soo, loao, 20.6, and 2011 modes respectively. As seen Figure [3] there are no global Alfven modes with 
frequency less than that of the 2ao mode. Thus, the values in Table U] correspond to the maximum magnetic field strength 
for each stellar model, if the observational data could be explained by the global Alfven modes. The maximum magnetic field 
strength is B/B^ 2 for SGR 1806-20 and B/B^^Z for SGR 1900-^14. 
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Figure 8. Comparison of several specific Alfvcn modes (solid lines) with observed frequencies for SGR 1900+14 (horizontal dashed lines) 
for the stellar model A+DH14. The agreement for all frequencies up to 150 Hz is excellent, for a magnetic field strength of -B/B^ « 1.94 
(vertical dashed-dotted line). 



5 SUMMARY AND DISCUSSION 

We have investigated torsional Alfven modes of relativistic stars with a global dipole magnetic field. It has been noted recently 
(Glampedakis et al. 2006) that such oscillation modes could serve as as an alternative explanation (in contrast to torsional 
crustal modes) for the SGR phenomenon, if the magnetic field is not confined to the crust. We compute global Alfven modes 
for a representative sample of equations of state and magnetar masses, in the ideal MHD approximation and ignoring ^ ± 2 
terms in the eigenfunction. The presence of a realistic crust has a negligible effect on Alfven modes for _B > 4 x 10^^ G. Using 
a toy model with large crust thickness we demonstrate that strong avoided crossings between torsional Alfven modes and 
torsional crust modes take place. For magnetar-like magnetic field strengths, the spacing between consecutive Alfven modes 
is of the same order as the gap of avoided crossings. As a result, it is not possible to identify modes of predominantly crustal 
character and all oscillations are predominantly Alfven-like. 

The Alfven-like character of global modes allows for a non-degeneracy of overtones with different i. Using this fact, we 
find excellent agreement between our computed frequencies and observed frequencies in two SGRs, for reasonable magnetic 
field strengths. We construct an empirical formula which describes the frequency of the fundamental torsional Alfven modes 
(for different €) as a function of the magnetic field strength and the compactness of the star. Such formulae could be used for 
iifferring the stellar properties and magnetic field strength. The requirement that the lowest observed frequency is identified 
with the fundamental 2^0 mode leads to a maximum magnetic field strength in the range of (0.8-1.2) x 10^® G. 

In the present work we have ignored couplings to ^ ± 2 terms. A two-dimensional treatment of the eigenvalue problem 
will appear elsewhere. 
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